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RESCUING THE INTRACLUSTER MEDIUM OF NGC 5813 


Noam Soker^, Shlomi HilleL, and Assaf Sternberg^ 

ABSTRACT 

We use recent X-ray observations of the intracluster medium (ICM) of the galaxy 
group NGC 5813 to confront theoretical studies of ICM thermal evolution with the 
newly derived ICM properties. We argue that the ICM of the cooling flow galaxy group 
NGC 5813 is more likely to be heated by mixing of post-shock jets’ gas residing in hot 
bubbles with the ICM, than by shocks or turbulent-heating. Shocks thermalize only 
a small fraction of their energy in the inner regions of the cooling flow; in order to 
adequately heat the inner part of the ICM, they would overheat the outer regions by a 
large factor, leading to its ejection from the group. Heating by mixing, that was found 
to be much more efficient than turbulent-heating and shocks-heating, hence, rescues the 
outer ICM of NGC 5813 from its predestined fate according to cooling flow feedback 
scenarios that are based on heating by shocks. 

Keywords: galaxies: active galaxies: clusters: general galaxies: groups: individual (NGC 
5813) - galaxies: jets 


1. INTRODUCTION 


A negative feedback mechanism driven by active galactic nucleus (AGN) jets that inflate X- 
ray bubbles (X-ray deficient cavities) determines the thermal evolution of the intra-clust er medium 
(ICM) in the inner region s of cooling flow (CF) clusters and groups of galaxies (e.g., iPfrommei 


2013 


Perucho et al.ll2014l l. In the feedback cycle the ICM feeds the AGN, and the AGN launches 


jets that heat the ICM. 

A consensus is emerging that the feeding of the AG N is through cold clumps, in what is termed 


cent vears have put the cold feedback mechanism on a verv solid ground (e.g.. 

Revaz et al. 

2008: 

Pope 

2009: Pizzolato & Soker 2010: Edee et al. 2 

010: Wilman et al. 

2011 
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2011; 
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^011: Caspar! et al.ll2012a b: 

McCourt et al.ll2012: Sharma et al. 

2012: Faraee et al. 

2012; 

Wash et al 

2014: Baneriee Sharma 

2014 

McNamara et al. 

2014; 

Voit &: Donahue 

2015: 
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Voit et al.ll2015: 

Li et al.ll2ni5: 

Prasad et al. 

2015 

Russell et al.ll2015: 

Tremblav et al.ll2015: 

Fogartv et al. 




The processes by which jets heat th e ICM, on the other h and, is still in dispute. B ased on ob¬ 
servations of the galaxy group NGC 581,1 iRandall et alJ ([20111) and iRandall et alJ (1201.51) argue that 


shock s heat the ICM; hereafter shocks-heating process. This was put into question bv ISoker et al 


(|2013lb and is further studied in the present paper. Base d on deep X-ray observation of the Perseus 
and Virgo cooling flow clusters, IZhuravleva et al.l (120141 1 argue that the main heating process is via 
dissipation of ICM turbulence; hereafter turbulent-heating process. The third heating process is the 
mixing of coolin g ICM with hot shock ed je ts’ gas that reside in the hot bubbles; hereafter mixing¬ 
heating process. ICilkis &: Sokerl (|2012l l and iHillel &: Soked (|2014l l conducted numerical simulations 
and showed that t he mixing-heat i ng pr ocess is much more efficient than the shocks-heating process. 
In a recent paper iHillel fc Sokerl (|201 (tI 1 showed that the mixing-heating is also more efficient than 
the turbulent-heating when turbulence is dr i ven by ACN activity. The turbulence, as observed 


in some clusters (e.g., IZhuravleva et al.ll201,5l : lAnderson fc Sunvaevll201,5h . is a consequence of the 


same vortices that cause mixing. 

To reach the conclusion that mixing-heating is the dominant process, it is mandatory to inflate 
bubbles by jets in a self-consistent manner. For example, the vortices that are induced inside and 


their evolution, and their i nteraction with the ICM, (e.g. 


2005; 


Roediger et al. 2007 : Sternberg fc Soker 2008b; Cilkis fc Sokerl 


(e.g. 

Omma et al. 

2004; 

Heinz &; Churazov 

Cilkis & Sokerl 

2012; 

Walg et al. 201, lib To 


inflate bubbles that resemble observed bubbles, our group has employed eith er slow (sub-relativistic) 


(Sternberg & Sokei 

2008a 

); see also 

observations (e.g.. 

Moe et al. 

2009: 


Falceta-Concalves et al. 


Dunn et al 


20101 : 


Arav et al 


e.g.. (ISternberg et al. 1 1200711 . or processing jets, e.g.. 


(I 2 OIO). SM W jets are supported by 


2013l b Both types of inflated 


bubbles will be described in the present study. 


In a recent paper iRandall et al.l (|201,5l j presented new observations of NCC 5813. Such detailed 
observations serve as a unique opportunity for us to confront our findings, regarding heating the 
ICM with jet-inflated bubbles, against observations. We here confront the new observations with 
two types of previously conducted numerical simulations. In section [3| we present a new study of 
the implications on the heating process of large scale flows that are induced by the jet-inflated 
bubbles. In section U] we show that a careful treatment is required when studying the shocks that 
are induced by jet-inflated bubbles. To our best knowledge, this is the first study that point to the 
infiuence of sound waves on the derived shock properties. But we first present in section [2] a new 
analysis of the effect on many repeated shocks on the outskirts of the ICM in NCC 5813. In section 
Owe summarize our findings, and conclude that although shocks play some role, mixing-heating is 
likely to play a much more significant role in the heating process. 
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IMPLICATIONS FOR THE OUTER ICM 


The inefficiency of the shocks-heating process was discussed in ISoker et alJ (j2013l ) . To excite 
a shock with a Mach number of Afg at radius r, the energy supplied by the AGN directly to the 
ICM, using the spherically symmetric Sedov solution for the shock, should be 

Fth(r), (1) 


where £'th(r) is the thermal energy of the ICM inside radius r. According to lRandall et alJ (|2ni5l '). 
to heat the cooling flow region of NGC 5813 inside 30 kpc, about 140 shocks with a Mach number 
of Aia ~ 1-2 are required during one cooling time. By equation ([1]), to reach this Mach number 
each shock is excited with an energy of Eghock ~ (l-2/l.l)^Ein = 1.2£'in, where we define = 
Eth{^ < 30 kpc). Each shock supplies an energy of Ein/140 = O.OlEin. Removing this dissipated 
energy in the inner region (up to 30 kpc), the energy that flows to the outer ICM regions within 
one cooling time at r = 30 kpc is Eout ~ 140(1.2 — 0.01)Ein ~ lOOEin. 

The ICM mass residing out at r > 30 kpc can be crudely estimated from the new density 
profile given bv iRandall et al.l (|2015l h The electron density profile from r = 22 kpc to r = 35 kpc, 
the largest radius in their profile, can be fitted with 

- 2.8 


ne(r) = 2.4 x 10 


-3 


cm 


-3 


( 2 ) 


30 kpc^ 

Extrapolating this fit to large distances of hundreds of kpc, we can estimate the mass in the outer 
region to be Mn„t ~ « 2Min. Based on ROSAT observations of poor groups of galaxies 


Mulchaev et al.l (| 19961 1 find the diffuse X-ray gas in a typical poor group to extend up to ~ 300 kpc, 
and have a shallower density profile. If we take a shallower density profile of 

- 1.5 


ne-s(r) = 2.4 x 10 


-3 


cm 


-3 


^30 kpc 

(like in eq. [5] below), the gas mass in the outer region to a distance of 300 kpc is Mout-, 
5 X IO^^Mq lOMin 


(3) 


Randall et al.l (|2015l l do mention that the energy carried by the bubble is transferred out by 


the buoyant bubbles. But the shocks turn into sound wave that dissipate their energy after a long 
time in the ICM. By the simple estimate made above the temperature of the outer gas should be 

(lit) 

This is a temperature larger than the virial temperature of the group. 

It seems that for the shocks-heating process to be effective in the inner region, the outer ICM 
would escape in a time much shorter than the cooling time, hence destroying the ICM observed 
structure. Clearly, a much more efficient heating mechanism than shocks is r equired to rescue the 
ICM of NGC 5813. Mixing-heating is this process ( Hillel fc Soker 2014 . 2016). 
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LARGE SCALE FLOW 


In the shocks-heating scenario presented by iRandall et alJ ([20111) and iRandall et alJ (j2015l ) 
tens of jet-launching episodes along the same axis heat the ICM. They do not consider any flow 
that might develop in the ICM. We here show that a large-scale meridional flow develops as a 
result of multi-jet activit y cycles. For that w e present results from 3D numerical simulations we 
have carried out recently ([Hillel fc Sokeiil2nR)l b We here briefly present the main characteristics of 
the numerical set up; more details can be found in that paper. 


We use the pluto code (jMignone et al.l 120071 1 for the hydrodynamic simulations in a three- 
dimensional Cartesian grid with adaptive mesh refinement (AMR). The computational grid is in 
the octant where the three coordinates x, y and z are positive. At the x = 0, y = 0 and z = 0 
planes we apply reflective boundary conditions. The z coordinate is chosen along the initial axis of 
the jets. The computational grid spans the cube 0 < x, y, z < 50 kpc. 

At the boundary z = 0 we inject into the grid a jet with a half-opening angle of = 70° 
([Sternberg et al. 1120071 1. The jet is injected during each active episode lasting 10 Myr, and when 
the jet is turned off for 10 Myr, reflective boundary conditions apply in the whole z = 0 plane. The 
initial jet velocity is uj = 8200 km s“^ corresponding to a Mach number of about 10 relative to the 
ICM. The power of the two jets together (we simulate only one jet) is E 2 j = 2 x 10^^ erg s“^ (half 
of it in each direction), and the mass deposition rate is M 2 j = 94M0 yr“^. 

The simulations presented here assume jets starting with wide opening angles. As stated in 
section [H there are observations of slow massive wide outflows from AGNs. A s well, the effect of 
rapidly processing jets in inflating bubbles is very similar to that of wide jets ([Sternberg &: Soker 
2008alb but it is simpler to simulate wide jets (for processing jets see next section). The main issue 
is to inflate bubbles, a process that involves the formation of vortices. Although magnetic fields 
are present in jets, as inferred from the radio emission, we here assume that they do not play a 
dynamical role. As we heat by mixing, we attribute no role to heat conduction on large scales (only 
on very small scales as mixing brings the ICM and the hot bubble gas very close to each other). 

The simulation begins with an isot hermal box of gas at an ini tial temperature of Ticm(O) = 
3 X 10^ K with a density profile of (e.g., IVernaleo &: Revnoldsll2006l ) 


/Oicm(?’) = 


Po 


1 -h (r/a) 


13/4 ’ 


( 5 ) 


with a = 100 kpc and po = 10“^^ g cm“^. A gravity field is added to maintain an initial hydrostatic 
equilibrium, and is kept constant in time. We include radiative cooling in the simu lations, where 
the tabulated cooling function is taken from Table 6 in [Sutherland fc Dopital (|l903l '). 


To reveal the large scale flow we follow an artificial flow quantity called ’tracer.’ The tracer 
is frozen-in to the flow, and hence tells us on the spreading with time of gas starting in a certain 
volume. A tracer’s initial value is set to .^(0) = 1 in a certain volume and ^(0) = 0 elsewhere. As the 
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traced gas mixes with the ICM or the jet’s material, its value drops to 0 < ^(t) < 1. We choose to 
trace the gas starting inside a torus around the z axis. The radius of the cross section of the torus is 
2.5 kpc, and it is centered at (y, z) = (20,15) kpc; or more generally + z) = (20,15) kpc. 


In Fig. [T]we present the evolution of the tracer at six times in the meridional plane (y, z). Note 
that the jets’ axis is along the z axis, which is the horizontal axis in the figure. T he evolution of 
other quantities in the (x, z) meridional plane is presented bv lHillel &: Soked (|201a i. In the panels 
on the left column the arrows represent flow velocity, with length linear with velocity up to an 
upper limit of Vm = 400 km s“^. Faster regions are presented with an arrow length corresponding 
to Vm- In the panels on the right column the arrows represent mass flux (p = vp. 


There are several prominent properties of the flow that are relevant to the present study and 
are revealed in Fig. [1] 


1. Mixing. The vortices created during the inflation of the bubble mix hot bubble gas with the 
ICM gas. This can be clearly seen by f ollowing the evolution of the tracer. This i s the main 


heating mechanism of the ICM by jets (ICilkis fc Soker 


2012; 


Hillel fc Sokeiil2014l . l2016l ). 


2. Outward then inward motion. The shocks running through the ICM followed by sound waves 
behind them can push gas out. The traced gas presented here is pushed from a torus cross 
section centered at (y, z) = (20,15) kpc to an elongated shape centered at (y, z) ~ (29,15) kpc 
at t ~ 100 Myr. However, then it turns around and flows i nward. It is heated o nly when it 
starts to be mixed with the hot bubble gas at t ~ 110 Myr (jHillel fc Sokeiil201(Tl ). 


3. Large scale meridional flow. At late times of t > 200 Myr, and after about 10 jet-activity 
episodes (cycles), a large-scale inflow takes place near the equatorial plane. This results 
from the dragging of ICM gas by the outflowing jets along their axis, here the horizontal 
z-axis. It seems though that in reality the jets’ axis will change direction and no such large 
meridional flo w will be develop. Th is change in direction can be seen in the three bubble-pairs 
in NGC 5813 ( Randall et al. 2015). 


It should be noted that the simulations presented here are general, and have parameters that 
fit clusters of galaxies, and do not fit groups of galaxies, such as NGC 5813. For example, the 
power of jets here i s 200-1000 times the power inferred for the bubble formation in NGC 5813 
(|Randall et al.ll201.5l b The consequence is that the mixing region extends much beyond a radius of 
30 kpc. For a much lower power of the jets as appropriate for NGC 5813, we would get a much 
smaller mixing region, one that fits the cooling flow region of NGC 5813. 


The main conclusion of this section is that one cannot assume that jets activity repeats itself 
along the same direction for many episodes, while the ICM in the equatorial plane (the plane 
perpendicular to the jets’ axis) does not flow inward. The flow inward will mix ICM gas with 
bubble gas, leading to heating by mixing that is much more efficient than the shocks-heating 
process. 
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Fig. 1.— Evolution with time of a traced gas simulated bv iHillel &: Soked (|2016l l. At t = 0 the 
traced gas was contained in a torus whose axis is the z-axis and whose cross section is a circle 
centered at z) = (20,15) kpc with a radius of r = 2.5 kpc. Color coding is of the tracer 

value The jets’ axis is along the horizontal z-axis in the panels. In the three panels on the left 
the arrow lengths are linear with velocity, up to a maximum value of Vm = 400 kms“^. Faster 
regions have an arrow length corresponding to Vm- The arrows in the right panels show the mass 
flux (p = vp. A length of 1 kpc on the map corresponds to (p = 2.3 x 10“^^ km s“^ g cm“^. 







































































































































































4. POST-SHOCK FLOW 


Randall et alJ (j2015l ) find the observed temperature rise across the shocks to be less than what 
is expected based on the Mach numbers they derive from their calculated density jumps across the 
shocks. They attribute the discrepancy to inability of the deprojected temperature measurements 
to resolve the tem perature jump due to the narrow width of the shock and the rarefied cool gas 
behind the shock (jRandall et al.l 1201 il l. 


We here raise anothe r possibility. When a bu bble is inflated by a jet sound waves are excited 
and trail the shock wave (jSternberg &: Sokeiil2009l ). The sound waves can compress gas behind the 
shock, and raise the density beh ind the shock much more than they raise the temperature. In Figs. 
[5] and E] we present results from ISternberg &: Soked ([20091), where all technical details are given. 


The relevant properties of these simulations are as follows. The simulati on are performed in 
half of the meridion al plane using the two-dimensional version of the VHl code (jBlondin et ahlllOOfll : 


Stevens et al.lll992l ) in spherical coordinates with a resolution of 256 x 256. Gravity was included 


in these simulations. Two opposite jets are launched along the x-axis. 


In the wide jets case two opposite jets were injected in the time period from t = —10 Myr 
until t = 0, each jet with a half opening angle of a = 70°. Their combined mass outflow rate was 
M 2 j = IOMq yr“^ and the initial velocity was radial with vj = 7750 km s“^. The total power of 
the two jets was E 2 ^ = 2 x 10^^ erg s“^. This case is presented in Fig. [2j 

In the case of the precessing jets the two opposite jets have the same mass outflow rate and 
velocity as in the wide-jets case, but they were injected at any given time with a half opening angle 
of a = 5°. In the axisymmetric grid employed, the jets are actually precessing very rapidly around 
the symmetry axis (i.e., in 3D torii were actually injected). The angle between the symmetry axis 
and the jets’ axis 9 was varied in a random way. The precession period, i.e., the time the jet returns 
to the same angle 6, is Tprec = 0.1 Myr. The jet’s interaction with the ICM is similar to that of 
a wide jet with a half opening angle of a ~ 50°. The jets were active between t = —18 Myr and 
t = 0. This case is presented in Fig. [3| 


The profiles along a cut in each case presented above show a density and temperature jump 
immediately behind the shock, the right edge in each of the two panels of profiles. Immediately 
behind the post-shock jump (moving to the left in the panels) there is another increase in the 
density, but the temperature decreases. The over all density jump overestimates the shock Mach 
number. 


We argue that iRandall et ahl (|2015l i somewhat overestimate the Mach numbers in the shocks 


they study in NGC 5813. 
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Fig. 2.— Excitation of sound waves by wide-jet-inflated bubbles. Upper panel: Density map for 
the ICM around the wide-jet-inflated bubbles at t = 10 Myr after the jets has ceased. The density 
units of the color-code bar are logp( g cm“^). The white regions are the two very-low density 
bubbles inflated by the two opposite jets. Lower panel: Density, temperature and pressure along 
the cut at 40° as marked on the upper panel by the dashed line. The pressure is given in units of 
6 X 10“^*^ erg cm“^, density in units of 1.5 x 10“^^ g cm“^, and temperature in units of 2.7 x 10^ K. 
Both panels are taken from Sternberg &: Soker (2009) where all details are given. The post-shock 
region at r = 21.5 kpc shows jumps in density, temperature, and pressure. The sound wave just 
behind the shock at r = 18.5 kpc has density and pressure above the post-shock values, while the 
temperature is lower than the post-shock value. 
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Fig. 3.— Excitation of sound waves by precessing-jet-inflated bubbles as simulated by Sternberg 
&: Soker (2009). Upper panel: Density map for the ICM around the precessing-jet-inflated bubbles 
at t = 0, just when the jets cease. Lower panel: Density, temperature and pressure along the cut 
at 80° as marked on the upper panel by the dashed-dotted line. The pressure is given in units of 
9 X 10“® erg cm“^, density in units of 2 x 10“^® g cm~^, and temperature in units of 2.7 x 10^ K. 
The post-shock region at r = 18 kpc shows jumps in density, temperature, and pressure. The sound 
wave behind the shock at r = 15.5 kpc has density and pressure above the post-shock values, while 
the temperature is lower than the post-shock value. 
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SUMMARY 


We reexamined the claim made by iRandall et alJ (j201ll ) and iRandall et alJ (j2015l i that the 
ICM in the cooling flow group of galaxies NGC 5813 is hea ted by many shock waves excited by 
repeated jet-launching episodes. It was already argued by ISoker et al.l (j2013l l that the shocks- 
heating mechanism is inefficient, and cannot account for the NGC 5813 ICM properties. We here, 
however, performed three new analyses and confronted theoretical studies with the newly derived 
properties of the ICM of NGC 5813. (i) We presented a new analytical analysis of the effect of 
many shocks on the outskirts of the ICM in NGC 5812. (ii) We discussed the role of the large 
scale flow induced by jet-inflated bubbles on the ICM and its thermal evolution, {in) We used 
simulations of sound waves trailing shocks induced by jet-inflated bubbles to deduce about possible 
complications in deriving the shocks’ properties from observations. These new analyses strengthen 
the conclusion that mixing-heating is more efficient that shock heating. 

The energy carried by the shocks and that is not deposited in the inner region, r < 30 kpc, is 
carried to the outer ICM regions and must dissipate there as sound waves. In section [2] we found 
that for the inefficient shocks-heating mechanism to heat the inner r < 30 kpc region, this extra 
energy heats the outer ICM to high enough temperature for it to escape the group of galaxies. If 
this was the case, the inner region would have expanded and escaped a long time ago. Clearly a 
more efficient heating mechanism is required. 

In section [3] we showed that after about ten jet-launching episodes along a constant axis a 
large-scale meridional flow develops (Fig{T]). An inflow of mass near the equatorial plane takes 
place. One cannot regard the ICM for such a case as static. Such a flow leads to a very vigorous 
mixing of the ICM with the hot bubbles’ gas. Heating by mixing becomes the dominate heating 
mech anism even for gas perpendicular to the jets’ axis (jCilkis &: Soked 120121 : |Hillel_&_Soker| [201J, 
2016l j. We expect that in reality the jets’ axis will change over time, and the meridional flow will 
not be as prominent as found here 

In section 0] we showed that sound waves propagating behind a shock wave can form a region 
behind the shock where the density in creases above the pos t-shock value, but the temperature does 
not. This might explain the finding of iRandall et al.l (|2ni5l i that the temperature jumps across the 
shocks in NGC 5813 are lower than those expected for the Mac h numbers they deduced based on 
the density jumps. If true, it implies that IRandall et al.l (|2015l l overestimated the Mach numbers 
in the three shocks they have found in NGC 5813. 

Our main conclusion is that shocks cannot be the main heating process of t he ICM in NGC 5813 . 
Heating by mixing s eems to be the main heating process of the ICM by jets (|Cilkis fc Sokeiil2ni2l : 
Hillel fc Sokerll2r)14l . 120161 1. The mixing is caused by vortices that are formed during the inflation of 
the bubbles. To reveal the properties of the vortices in numerical simulations it is crucial to inflate 
bubbles in a self-consistent manner. When that is done, it is found that the mixing operates also 
perpendicular to the jets’ axis. 
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The vortices lead to turbule nce in the ICM. The refor e, the finding of mod erate turbulence in 
the inne r region of NGC 5813 bvlWerner et alJ (1200911 andlde Plaa et al.l (j2012l b as well as in other 
clusters (jZhuravleva et alJ l2014l . l2015l : lAnderson &: Sunvaevi [201a) , strengthens the argument that 
the mixing-heating process is the dominate heating process. Some heating by the dissipation of 
turbu lence takes place as well, but it accounts for at most ~ 20% of the total heating (jHillel fc Soker 


201 Ol i. 
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